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ABSTRACT 

We present a new perspective on the crystallinity of dust in protoplanetary disks. The dominant 
crystallization by thermal annealing happens in the very early phases of disk formation and evolution. 
Both the disk properties and the level of crystallinity are thereby directly linked to the properties 
of the molecular cloud core from which the star+disk system was formed. We show that, under the 
assumption of single star formation, rapidly rotating clouds produce disks which, after the main infall 
phase (i.e. in the optically revealed class II phase), are rather massive and have a high accretion rate 
but low crystallinity. Slowly rotating clouds, on the other hand, produce less massive disks with lower 
accretion rate, but high levels of crystallinity. Cloud fragmentation and the formation of multiple 
stars complicates the problem and necessitates further study. The underlying physics of the model 
is insufficiently understood to provide the precise relationship between crystallinity, disk mass and 
accretion rate. But the fact that with 'standard' input physics the model produces disks which, in 
comparison to observations, appear to have either too high levels of crystallinity or too high disk 
masses, demonstrates that the comparison of these models to observations can place strong contraints 
on the disk physics. The question to ask is not why some sources are so crystalline, but why some 
other sources have such a low level of crystallinity. 

Subject headings: accretion disks — (ISM:) dust — (stars:) planetary systems: formation, protoplan- 
etary disks 



1. INTRODUCTION 

One of the major discoveries achieved with infrared 
spectroscopic studies of Herbig Ae/Be stars and classi- 
cal T Tauri stars is that a significant percentage of the 
dust in their circum stellar disks is of cryst alline form 
(e.g. Bouw man et al. 120011 van Boekel et al. 2005; Apai 
et al. I2005D . Also in our own solar system there is evi- 
dence, from observations of comets, that the dust of the 
primordial solar nebula was partly crystalline (Wooden 
et al. 0,999). This poses an interesting puzzle because the 
dust inherited from the interste llar medium must have 
been amorphous (Kemper et al. 2004; 2005) and hence 
the crystallization must have taken place within the disk. 
Since silicate dust crystallizes only when heated to rel- 
atively high temperatures (> 800 K), the discovery of 
these crystals proved that the dust in such disks is (or 
has been) subject to processes involving high tempera- 
tures. 

The very inner regions of protoplanetary disks (i.e. in- 
ward of the 'crystallization radius' which is about 0.1 
AU for a T Tauri star and 0.7 AU for a Herbig star) are 
warm enough to thermally anneal the dust, and thus pro- 
duce crystalline silicates like enstatite and forsterite. But 
the presence of crystalline dust in Solar System comets 
proves that processed dust also existed beyond the snow 
line where comets formed, well ou tside the c rystall ization 
radius (e.g. Wooden et al. 119991: Min et al. 120051) . Sim- 
ilarly, the infrared spectra of T Tauri and Herbig disks 



show that crystalline silicates often exist out to radii as 
cool as 150 K. And by directly spatially resolving the disk 
down to 2 AU with interferometric 10 /im spectroscopy 
it has been shown that for many sources the outer cooler 
disk region s clearly contain crystalline dust (van Boekel 
et al. I2004|) . though a smaller fraction than in the inner 
disk regions. 

Curren tly the most favored theory involves radial mix- 
ing ( Gail 120011 120021 Wehrstedt & Gail l2002t Bockelee- 
Morv an et al. 120021 based on earlier ideas by Morrill & 
Voelk ll984H . Even though the main accretion stream in 
a disk points inward toward the star, the turbulent mix- 
ing within the disk can transfer at least a small amount 
of the thermally processed dust grains to larger radii. 
This typically results in decrease of crystalline silicate 
abundance with radius beyond the point where the tem- 
perature drops below the annealing temperature. How- 
ever, this decline is slow enough that still a measure- 
able quantity of crystalline silicates exists at radii where 
the 10 /im silicate feature is produced in protoplanetary 
disks, and where comets were f ormed in the protosolar 
nebula (Bockelee-Morvan et al. I2002D . Whether or not 
thermal annealing and radial mixing alone are sufficient 
to explain the levels of crystallinity typically observed is 
debated. A possibility to enhance the outward transport 
of cry stallin e silicates was recently discussed by Keller & 
Gail (2004), who have shown that accretion in the in- 
ner disk regions is not always inward: While the surface 
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layers of a disk tend to accrete inward, the midplane of 
the disk moves outward. Any thermally processed grains 
may be transported very efficiently outward with this 
mechanism. Another possibility is t he ou tward trans- 
portation of dust via wind (Shu et al. Il99fith 

In this Letter we wish to revisit the radial mixing the- 
ory from a new perspective. Instead of assuming a steady 
disk, or starting from a given disk structure, we take a 
step back and start from the collapse of the rotating pre- 
stellar cloud core. We assume that the collapse of such 
a core leads to the formation of a single star and a disk, 
even though this may not be the dominant mode of star 
formation. The properties and viscous evolution of this 
disk are strongly influenced by the rotatio n rate of the 
original cl oud co re (Nakamoto & Nakagawa Hueso 
& Guillot 2005). Rapidly rotating clouds deposit their 
mass onto the disk at large stellocentric radii, leading 
to the slower evolution of a more massive disk than if 
the cloud were rotating slowly. Slowly rotating clouds 
deposit most material close to the star where it will 
be strongly heated and then push outwards by viscous 
spreading, influencing the thermal history of material 
found in th e outer disk (Fig. 0J. Naka moto & Naka- 
gawa ()1994fl and Hueso & Guillot ( 2005) have modeled 
the formation and viscous evolution of protoplanetary 
disks using vertically integrated time-dependent viscous 
disk models. Here we follow a similar approach but along 
with the disk evolution we also passively evolve a dust 
population, including a simple treatment of crystalliza- 
tion and radial mixing. 

In this Letter we only demonstrate the principle and 
present first results. In a follow-up paper we will describe 
the model in detail, explore the parameter space, discuss 
all the physical uncertainties and caveats, and compare 
the results to observational trends. 

2. BRIEF DESCRIPTION OF THE MODEL 

We start our simulations with a core with tempera- 
ture Tbg (Kelvin), mass M (gram), and solid-body ro- 
tation rate O (radian/second). For the sake of sim- 
plicity we assume the core to be a singular isothermal 
sphere, required to apply the Shu collapse model. Here 
the collapse is triggered by an expansion wave propagat- 
ing outward from the center with the sound speed c s = 
\/kT\ }% j \im v , initiating the collapse of every mass shell 
of the core it passes through. The infall rate is constant: 
Onfall = 0.975 c^/G, which for our choice of tempera- 
ture, T bg = 15 K, amounts to Mi n f a ii = 3x lO~ 6 M /year. 
For a cloud with a rotation rate fl all the mass will fall 
onto the equatorial plane inward of the centrifugal radius 
i? contr = il 2 R^ OIC /GM COIC , where the radius of the core 
is i? C orc — GAI colc /2c 2 . The way the infalling material is 
distrib uted over the disk is described by Hueso & Guillot 
(20QU). 

The disk is described by the surface density X as a 
function of the radius R. It obeys the equation: 

as 1 dKZvR _ 

~m + R~^R~~ b ' (1) 

where S is the source function due to infalling matter 
onto the disk. The azimuthal momentum equation yields 
an expression for the radial velocity of the gas vr: 
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Fig. 1. — Surface density as a function of radius in the disk, for 
different times after the onset of collapse of the parent cloud. 

where v is the viscosity coefficient, given by v = 
akT m / '[j,m p flK- He re a is the parameter of viscosity 
(Shakura & Sunvaev ll97.il) . which we take to be a global 
parameter, T m is the midplane temperature, and £Ik is 
the Keplerian frequency. We also include an effective 
viscosity caused by gravitational instabilities in regions 
where the Toomre parameter drops below unity. The 
midplane temperature is determined by taking into ac- 
count heating due to viscosity and irradiation by the cen- 
tral star. 

The dust is passively transported along with the gas, 
but is subject to radial diffusion through turbulence and 
to crystallization in those regions for which T m > 800 K. 
We therefore transport two dust components: an amor- 
phous one and a crystalline one, with surface densities 
Si and £2, respectively. Crystallization transforms one 
into the other at a rate given by the Arrhenius formula 
1/ exp(—T c /T m ) s _1 with t he latt ice vibration frequency 
v = 2 x 10 13 (Fabian et al. l2005l an d T c = 38100 K (ta- 
ble 3 of Bockelee-Morvan et al. 2002). The radial mixing 
equation is: 

dZi 1 dREiV R _ 1 d 
~dT + R dR ~ RdR 



RDE± ( * 
dR V £ 



Si (3) 



where D = v/P t is the diffusion coefficient, Pt is a con- 
stant that follows from numerical simulations (Johansen 
& Klahr l2005l Carballido et al. l2005jl . The source term 
Si accounts for the infall of new matter onto the disk 
(only feeding the amorphous component) and for the 
transformation of amorphous silicates into crystalline 
ones. We assume that the amorphous primordial dust 
survives the accretion shock on the surface of the disk. 

3. RESULTS 

As an example, we take a Herbig Ae/Be star with 
M* = 2.5M Q , T* = 10,000 K, L = 50 L Q . The ini- 
tial cloud is assumed to be at 15 K, and its rotation rate 
is taken to be Q = 1 x 10 _14 s _1 . We assume that the 
viscosity parameter is a = 0.01 throughout the disk. 

In Fig. ^ the surface density profile as a function of 
radius is shown at various epochs after the onset of col- 
lapse of the parent cloud. In this example the centrifugal 
radius was located at 180 AU, i.e. all the matter from 
the infalling cloud fell within 180 AU from the star onto 
the disk. Viscous friction then causes much of the mass 
to accrete inward onto the star, while pushing a smaller 
fraction of the matter out to very large radii in order 
to absorb the angular momentum of the inward moving 
matter. 
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Fig. 2. — Abundance of crystalline silicates as a function of 
radius, for the same times as in Fig. HI 

In the initial phases of the infall (t <C 10 6 year), how- 
ever, the centrifugal radius of the infalling matter was 
smaller than 180 AU because the inner parts of the 
rigidly rotating cloud had a lower specific angular mo- 
mentum. Therefore, in the early phases of the collapse 
most of the matter falls much closer to the star than in 
the later stages. In these inner regions the disk becomes 
so warm due to accretion that all the dust crystallizes 
or evaporates. As some of this matter is pushed out for 
reasons of angular momentum conservations (not nec- 
essarily radial mixing), this matter cools off and, if the 
dust was evaporated initially, it recondenses directly into 
crystalline form. Hence, the spreading disk is transfer- 
ring crystalline matter from the inner disk to large radii, 
causing initially the disk to be fully crystalline, as shown 
in Fig. [21 

In the later stages of the infall (~ 5x 10 5 • • ■ 10 6 years), 
the collapse wave reaches the outer shells of the cloud, 
triggering this to fall onto the disk. Since this matter 
has a higher specific angular momentum it falls onto 
the disk at larger radii, where the disk is cooler. This 
matter therefore will stay amorphous and mix with the 
crystalline material already in the disk, diluting the 
crystalline content. The disk becomes more amorphous 
again. The onset of this dilution with amorphous ma- 
terial is seen in the 3.3 x 10 5 year curve (dotted curve) 
in Fig. |3 Dilution continues until infall ceases, which is 
at slightly less than 1 Myr in this model. Further vis- 
cous spreading and radial mixing of the dust will cause 
the region beyond 10 AU to have an almost constant 
crystallinity, while in between 1 and 10 AU crystallinity 
declines with radius, in accordance with the radial mix- 
ing theory. The level of the crystallinity in the outer 
(10-1000 AU) disk region is a heritage from the infall 
phase set by the amout of cool non-crystallized matter 
that was able to dilute the crystalline matter from the 
very early infall phase. 

In Fig. 01 we plot the crystallinity at 10 AU at 2 Myr 
after the onset of collapse for a series of models with 
different cloud rotation rate Vl. Also shown are the disk 
masses and accretion rates at 2 Myr. 

The figure shows that cloud rotation rates of at least 
VL > 4 x 10~ 15 s~ 1 are required to obtain crystallinity 
below 30%, which is the approximate upper limit ob- 
tained for a s ample of Herbig Ae/Be stars by van Boekel 
et al. (|2005l) . Disks formed from such rapidly rotat- 
ing clouds are large and require long time to be ac- 
creted. At 2 Myr they will still have high disk masses 
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Fig. 3. — A number of disk properties plotted as a function of 
cloud rotation rate f2, or equivalently as a function of the centrifu- 
gal radius i? C cntr- The latter is the fundamental variable relevant 
for the disk evolution, but with the simple cloud collapse model 
it translates in O. The following quantities are plotted: the crys- 
tallinity at 10 AU (i.e. abundance of crystalline silicates: black 
solid curves), the mass accretion rate at the inner edge of the disk 
(dotted curve), and the disk mass (dashed curve), all measured at 
2 Myr after onset of collapse, for the model parameters described 
in the text. The top black curve is for the model with the "stan- 
dard" input physics. The bottom black curve is for the model in 
which the most optimistically low midplane temperature estimate 
is used, and in which the radial mixing efficiency is reduced by a 
factor of 10 (see text). 

(A'idisk ^ 0.4M Q ) and reasonably high accretion rate 
(M > 10" 7 M^/year). 

The crystallinity depends strongly on the detailed 
treatment of the temperature at the disk midplane. Our 
standard computation leads to rather hot disks, but 
many effects might lower the disk temperature, includ- 
ing convection and dust grain growth. The minimum 
temperature is the effective temperature required to ra- 
diate away the energy input from both the irradiation 
as well as the viscous heating. This lower temperature 
will evidently yield lower crystallinity. Crystallinity can 
be suppressed even further by reducing the radial mix- 
ing coefficient (i.e. increasing the Prandtl number P t ). 
The grey curve in Fig. represents models with P t — 10 
and minimum temperature. For these models a level 
of 30% crystallinity can be reached (Mdi s k =0.05 M Q , 
M = lO~ 8 M /year at 2 Myr). 

However, the above assumptions differ drastically from 
conventional wisdom: The lowest possible estimate of 
the disk midplane temperature would require a Rosse- 
land optical depth of roughly unity throughout the disk, 
a rather unlikely coincidence. The real temperature pre- 
sumably lies somewhere in between the two extremes we 
tested. The factor of 10 reduction in radial mixing means 
that the angular momentum in the disk is much more ef- 
ficiently transported than passive trace elements in the 
disk. In accretion disk theory the Prandtl number is typ- 
ically assumed to be about unity. 

4. DISCUSSION AND CONCLUSION 

In this Letter we have shown that understanding the 
crystallinity of dust in disks requires the modeling of the 
entire formation process of the star+disk system, be- 
cause much of the crystallization occurs in these very 
early phases. Unlike the existing models in the literature 
where the initially amorphous disks are gradually en- 
riched with crystals via radial mixing, our model predicts 
a major heating event and dust crystallization simulta- 
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Fig. 4. — Cartoon of how the angular momentum of the infalling matter affects the crystallinity of the resulting disk. Left: low angular 
momentum infalling matter falls onto the disk within a small region around the star. In this region the disk (in particular in the early active 
phase) is hot enough to crystallize the du st. Subsequent disk-expansion transports this material outward (some of which later accretes 
back inward, Lynden-Bell 8z Pringle Il974l) . Much (if not most) of the matter observed at large radii originally came from these hot inner 
regions, and hence has a high abundance of crystalline silicates. Right: high angular momentum matter falls onto the disk at large radii, 
where the disk is cool. This matter will not get crystallized. Even though radial mixing will pollute it partially with crystalline silicates 
from smaller radii, the overall crystallinity of the disk will then be much lower than in the high-angular-momentum case. 



neously with the build-up of the accretion disk. Thus, 
in contrast to the gradual increase of crystallinity pro- 
posed by others, our model predicts an initially high crys- 
tallinity, that possibly may even decrease with time — a 
prediction that might provide a natural explanation for 
the probable age-cry st allin it y anti-correlat ion in dicated 
by van Boekel et al. l)2005j) and Apai et al. 120051 

The crystallinity, as well as many other disk parame- 
ters measured after the main envelope-infall phase, de- 
pends crucially on the angular momentum of the in- 
falling matter: low angular momentum leads to high 
crystallinity and high angular momentum to low crys- 
tallinity. Our model only includes thermal annealing as 
a crystallization process. Other processes such as shock 
annealing (both within the disk itself and in the infall- 
shock on the disk surface) are neglected. Nevertheless, 
our model typically predicts high levels of crystallinity. A 
detailed comparison to observations (e.g. to Spitzcr IRS 
spectra and/or VLT-MIDI interferometric spectroscopy 
data) can therefore put strong constraints on the disk 
physics. The low levels of crystallinity observed in some 
sources may require lowering the vertical optical depth 



of the disk by dust coagulation. We may also be forced 
to reduce the radial mixing or modify other pieces of 
basic disk physics. The study of the disk crystallinity, 
in combination with disk mass and accretion rate, will 
therefore teach us about several aspects of fundamental 
disk physics. In upcoming work we will study these as- 
pects in more detail and perform direct comparisons to 
observations. 
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